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Abstract 

A solution to the BBGKY hierarchy for nonequilibrium distribution functions is obtained within 
modified boundary conditions. The boundary conditions take into account explicitly both the non- 
equilibrium one-particle distribution function as well as local conservation laws. As a result, modified 
group expansions are proposed. On the basis of these expansions, a generalized kinetic equation for 
hard spheres and a generalized Bogolubov-Lenard-Balescu kinetic equation for a dense electron gas 
are derived within the polarization approximation. 
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1 Introduction 

The construction of kinetic equations for dense gases, liquids and plasma, where a small parameter 
(density or interaction) is absent, still remains one of the most significant complication in the kinetic 
theory of classical systems. A number of difficult tasks centre around it. The study of solutions to the 
BBGKY hierarchy for nonequilibrium distribution functions in view of interparticle correlations could be 
considered as a way to solve this problem. 

Recently, a kinetic equation of the revised Enskog theory for a dense system of hard spheres and an 
Enskog-Landau kinetic equation for a dense system of charged hard spheres have been obtained from 
the BBGKY hierarchy in the "pair" collisions approximation jlL 0, t|. It is necessary to notice that this 
approximation does not correspond to the usual two-particle approximation, inherent in the Boltzmann 
theory, because an essential part of the many-particle correlations is implicitly taken into account by the 
pair quasiequilibrium distribution function g2{ri, T"2;i). 

To analyze solutions to the BBGKY hierarchy ^ in higher approximations on interparticle correla- 
tions, it is more convenient to use the concept of group expansion s 0, |q, pL This has been applied to the 
BBGKY hierarchy in previous investigations j§ |, |, |, |, |, g ]13|, [[J, |l§ [[§ [lT| using bound- 

ary conditions which correspond to the weakening correlations principle by Bogolubov |18|. The same 
conception has been involved in papers of Zubarev and Novikov [l^, |2^, |2l[ , where a diagram method 
for obtaining solutions to the BBGKY hierarchy was developed. Moreover, in papers by D.N. Zubarev et 
al. [^2] , p3|] a consistent description of kinetics and hydrodynamics has been proposed. It is based on a 
new formulation of boundary conditions for the Liouville equation and BBGKY hierarchy. In particular, 
this formulation takes into account interparticle correlations which are connected with local conservation 
laws. 

In this paper, the conception of group expansions will be applied to the BBGKY hierarchy with 
modified boundary conditions which take into account both the nonequilibrium behaviour of one-particle 
distribution function and local conservation laws (sec. 2). In sections 3 and 4 it will be demonstrated on 
the basis of the modified group expansions how to obtain a generalization of some known kinetic equations 
for hard spheres and homogeneous plasma to the case of high densities. Advantages and shortcomings of 
the proposed approach are discussed concisely in conclusion and at the end of some sections. 



2 The BBGKY hierarchy with modified boundary conditions 



The BBGKY hierarchy of equations for nonequilibrium distribution functions of classical interacting 
particles has been obtained in the paper |2^] on the basis of assembling of time retarded solutions for the 
Liouville equation with the modified Bogolubov's condition of weakening correlations between particles. 
According to the Zubarev's nonequilibrium statistical operator method |24|, ^5fl , the total nonequilibrium 
distribution function g (a; ; t) for all N particles of a system satisfies the following asymptotic condition: 



lim exp(iL N t )(g(x N ,t ) - g q (x N ,t )) = 0. (2.1) 



Here, i = \J — 1, and the limit to — > — oo is taken after thermodynamic limiting transition N 
V — > oo, N/V — > const; Lm is the Liouville operator: 

j=i j=i fe=i ' v * • ' •• 



is the interaction energy between two particles j and fc; = is the set of phase variables 

(coordinates and momenta). 

Let us consider time scales when details of initial state g [x N , to) become not so important, i.e., when 
t to. Then, to avoid the dependency on to, let us average the formal solution to the Liouville equation 
g(x N ;t} = c -' L ( t -*o)g (x N ; tp) with respect to initial times in the range to ^ t and make the limiting 
transition to — t — > oo [E4| E5[ , 



g(x N ;t)=e / dt' e £t 'e iLNt 'g q (/; t + t') , (2.3) 



where e tends to +0 after the thermodynamic limiting transition. The quasiequilibrium distribution 
function g q (x N ; t) is determined from the condition of maximum for informational entropy at fixed values 
of the single particle distribution function fx (xx; t) and average density of the interaction energy (£i n t(r))', 
(...)* = / dPiv • • • Q (x N ;t), d_Tjv = (dx) N /Nl. This corresponds to taking into account correlations, 
related to the conservations laws of hydrodynamical variables for particles density n(r;t), momentum 
j(r; t) and full energy £(r;t) 



The boundary condition (2.1) is equivalent to the transition from the Liouville equation to a modified 
one @|5|: 

(j t + iL N ^j g (x N ,t) = -s(g (x N ,t) g q (x N , t) ) . (2.4) 

This equation contains a small source on the right-hand side, which destroys the invariance with respect 
to time inversion. 

The choice of g q (x N ; t) depends mainly on nonequilibrium state of the system under consideration. In 
the case of low density gases, where times of free motion are essentially greater then collision times, higher- 
order distribution functions of particles become depend on time via oneparticle distribution functions 
only ||. It does mean that an abbreviated description of nonequilibrium states is available, and the total 
nonequilibrium distribution function depends on time via fx(x;t). In such a case, the quasiequilibrium 
distribution function g q (x N ; tj reads || : 

*(^;*)=n^^, (2.5) 



where e is the natural logarithm base. Then, the Liouville equation with a small source ( |2.4| ) in view 
of ( |2.5| ) corresponds to the abbreviated description of time evolution of the system on a kinetic stage 
when only the oneparticle distribution function is considered as a slow variable. However, there are 
always additional quantities, which vary in time slowly, because they are locally conserved. In the case 
of an onecomponent system, the mass density p(r; t), momentum j(r\ t) and total energy £(r; t) belong 
to those quantities. At long times they satisfy the generalized hydrodynamics equations. Generally 
speaking, the equation for fx(x;t) should be conjugated with these equations. For low density gases, 
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such a conjugation can be made, in principle, with arbitrary precision in each order over density. In 
high density gases and liquids, when a small parameter is absent, the correlation times corresponding to 
hydrodynamic quantities become to be commensurable with characteristic times of varying of oneparticle 
distribution functions. Therefore, in dense gases and liquids the kinetics and hydrodynamics are closely 
connected between themselves and they should be considered simultaneously. That is why, manyparticle 
correlations, related to the local conservation laws of mass, momentum and total energy, can not be 



neglected |2J, 25 . The local conservation laws affect on kinetic processes, because of the interaction 
of selected particles group with other particles of a system. This interaction is especially important 
in the case of high densities, and it must be included into consideration. Then at the construction of 
kinetic equations for high densities, it is necessary to choose the abbreviated description of nonequilibrium 
system in such a form to satisfy the true dynamics of conserved quantities automatically. To this end, 
the densities of hydrodynamic variables should be included together with the oneparticle distribution 
function fi(x;t) into the set of parameters of the abbreviated description initially |l], ^2], |23|. The next 
phase functions correspond to densities of hydrodynamic variables p(r; t), j(r; t) and £(r; t): 

p(r) = fdpmhi(x), £ (r) = /dp f^ni(x) + 1/2/dr' dp dp' ${\r - r'\)n 2 (x, x'), 
j(r) = fdpphi(x). 

where fii(x) and ^(x, x') are the one- and two-particle microscopic phase densities by Klimontovich JT3| |, 
Relations ( |2.6| ) show a distinctive role of potential interaction energy. In contrary to p(r;t) — (^(r))' 
and j(r;t) — (j(t"))*, nonequilibrium values of the total energy £(r;t) — (£ (r))* can not be expressed 
via oneparticle distribution function f\{x\t) = {h^x)) 1 exclusively, because to evaluate a potential part 
of £; nt (r;i) = (fi nt (r))* it is necessary to involve the twoparticle distribution function f 2 (x,x';t) — 
(n 2 (x, x')) 1 . Here 

&rt(r) - \ J dr' dp dp' $(|r - r'\)h 2 {x, x') (2.7) 

is the density of potential energy of interaction. The next conclusion can be formulated as follows. If 
oneparticle distribution function fi(x ; t) is chosen as a parameter of the abbreviated description, then 



the density of interaction energy (2.7) can be considered as an additional independent parameter. One 
can find the quasiequilibrium distribution function g q (x N ;tj from the condition of extremum for the 
functional of informational entropy Si n {(t) = — J dJV Q (x ; t) In g (x N ; t) at fixed average values of 
{hi{x)Y = fi(x;t), (£int( r )) t = £int( r ;*) including the normalization condition for g q (x N ;t). After 
simple transformations the desired relation for the quasiequilibrium distribution function reads 

g q (x N ;t) = exp|-$(t)- J dr f3(r;t)£ int (r) - J dx a(z; (a) | , (2.8) 

$(0 = lnjdr N expj- J dr (3(r;t)£ int (r) - J dx aix-^fnix)^, (2.9) 
where $(t), /3(r;t), a(x;t) are the Lagrange multipliers. <!>(£) is the Massieu-Planck functional, which is 



determined from the condition of normalization for g q (x N ;ty. J d/V^q (x ji) — 1- Relation (2.8) for 



g q (x N ; t) has been obtained for the first time in (2^] . To determine the physical meaning of parameters 



(3(r;t) and a(x;t) let us rewrite g q (x N ;t) (2_^8) in the form 

e , (*»;«) = =p{-*( f ) - / *. Rr-tfir) - / d, A^)} , (2.10) 

where £'(r) is the density of total energy in a reference frame, which runs together a system element 
with the mass velocity V(r;t) [^3[ Q 

£'{r) =£{r)- V(r; t)j(r) + - V 2 (r; t)n(r). (2.11) 



Here n(r) = J dp fii(x) is the density of particles number. Parameters f3(r;t) and a'(x;t) in (2.10) are 
defined from conditions of self-consistency, namely, the equality of quasiaveraging values (hi(x)) q and 
(£' (r))* to their real average values {rii(x)} t , (£'{r)) t : 

(niCc))« = (Wx))* = h(x;t), 
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Here (. . .)* = J dl^r • • • Qq (x N t) ■ In these transformations parameters a'(x; t) and a(x; t) are connected 
by the relation 



a'(x; t) = a(x; t) - /3(r; t) ( f- - V(r; t)p + ^V 2 (r; t)\ . 

[ 2m 2 J 



In the case, when conditions ( 2.12 ) take place, one can obtain some relations, taking into account the 
self-consistency conditions and varying the modified Massieu-Planck functional 



= In / dTjv exp 



dr f3(r;t)£'(r) - / dx a'(x;t)n 1 (x) 



with respect to parameters (3(r;t) and a'(x;t) 
6${t) 



sp(r,t) 

5a(x; t) 



= -<£'(r))*=-<£'(r)>*, 
= -{nx{x))\ \=-{n x {x)Y = -h{x;t). 



(2.13) 



It means that the parameter f3(r;t) is conjugated to the average energy in an accompanying reference 
frame, a'(x; t) is conjugated to the nonequilibrium oncparticle distribution function fi(x; t). To determine 
the physical meaning of these parameters let us define entropy of the system taking into account the self- 
consistency conditions ( p. 12 ) 



S(t) = -fln&j {x N ;t))l = $(t) + / dr /3(r; t)(£ '(r))' + / dx a'(x; t)<ni(o;)) 



(2.14) 



Taking functional derivatives of S'(t) ( ^.14| ) with respect to (£'(r))' and (ni(x))' at fixed corresponding 
averaged values gives the following thermodynamic relations: 



8(£'{r)Y oh(x;t) 



(2.15) 



Hence, (3(r;t) is an analogue of local inverse temperature. 

In general case, when kinetic and hydrodynamic processes are considered simultaneously, the quasi- 
equilibrium distribution function ( 2.10| ) or (|2.8[ ) can be rewritten in somewhat other form. This form is 
more convenient for comparison with g q (x N ; (2.5), obtained in the usual way |^|, when fi(x;t) is only 
the parameter of the abbreviated description. First of all, let us note that one can include the parameter 
$(t) in ( |2.8[ ) into parameter a(x; t) as a term, which does not depend on x. Parameter a{x; t) in g q [x N ; t) 



can be excluded with the help of the self-consistency condition (n(x)) q = (h(x)) t — fi{x; t). Reducing of 
g q (x N ; t) results in 



A' 



,.N 



t) =cxp{-U N (r N ;t)}]J 



i=i 



u(rr,t) ' 



where functions u(ri\t) are obtained from the relations 
u(ri;t) = 



drAr ^exp{-^(r,r-^)}n n(r ' :/) 



(N-iy. 



1 = 2 



u(ri;t) ' 



(2.16) 



(2.17) 



1 N 

U N (r N ;t) = U N ( ri ,...,r N ;t) = - ^ $(\ ri - r^r^t)., 



1^3 = 1 



n(r;t) — (ri.(r))* = J dp fi(x;t) is the nonequilibrium particles concentration. In last expression ( 2.16 ). 
/7jv(r;t) and u(ri;t) depend explicitly and implicitly, respectively, on n(r;t) and /3(r;t)(or (£'(r))*). To 
obtain the ordinary Bogolubov scheme H, it is necessary to put U~N(r;t) = in ( |2.16|) a nd (2.17). Then, 
one can define u — e, and (2.16) transforms into the quasicquilibrium distribution (2.5) as it should be. 
In general case u(r;t) is a functional of nonequilibrium density of particles number (n(r)) and (3(r;t), 
which is an analogue of local inverse temperature. Nevertheless, one should handle with care with this 
analogy as far as definition ( 2.16| ) can describe states which are far from local equilibrium. In particular, 
fi(x;t) can distinguish considerably from the local Maxwellian distribution. 
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The entropy expression (2.14) can be transformed accordingly to the structure of the quasiequilibrium 
distribution function (2.16) 



S(t) 



dr /3(r;t)(S in t(r)y 



dx fi{x; t) In 



u(r;t) ' 



(2.18) 



Here "potential" and "kinetic" parts are separated. In the case of low densi ty ga ses, the influence of 
the potential energy can be neglected and u(r;t) — e. Then the expression ( |2.1S| ) tends to the usual 
Boltzmann entropy. 



3 The modified group expansions 



Integrating equation (2.4) over the phase space of (N — s) particles, we obtain the equation chain for the 
s-particle nonequilibrium distribution function f a (x s ;t) — j d.Tjv- s f) (x N ;t^ [jl], p2[ : 



{jt + [Ls ) fs {xS;t) + ^ / dXs+1 [L{J ' S + 1)/s+1 



(3.1) 



-e y s (x s ;t)-g s (r s ;t) J[ h (xj 5 1) J , 

where g s (r s ;t) = f dJjv— s dp s g q (x N ,t) is the quasiequilibrium s-particle coordinate distribution func- 
tion, which depends on n(r;t) and /3(r;t) functionally. Due to the fact, that gi(ri;t) = 1, the equation 



chain (3.1) is distinguished from the ordinary BBGKY hierarchy |L8| by the existence of sources in the 
right-hand parts of the equations beginning from the second one. Such sources take into account both 
one-particle and collective hydrodynamical effects. 



Let us note, that equation chain (3.1) should be amplified by equations for spatial quasiequilibrium 
distribution functions g s (r s ;t), which functionally depend on nonequilibrium density of particles num- 
ber n(r;t) and "inversed" local temperature (3{r;t). Specifically, it has been shown in ^6|, that pair 
quasiequilibrium distribution function g2(r\,T%\t) is connected with a pair quasiequilibrium correlation 
function: h 2 (ri,r2;t) = 52(^1, ^2; t) — 1. On its turn, ^2(^1)^2;*) satisfies the Ornstein-Zernike equation 



h 2 {r 1 ,r 2 ; t) = c 2 (n, r 2 ; t) + J dr 3 c 2 (r 1 , r 3 ; t)h 2 {r 2 , r 3 ; t)n(r 3 ; t), 



where c 2 (r\,r 2 ] t) is a direct quasiequilibrium correlation function. 

To analyze the BBGKY hierarchy, we follow similarly to the papers by Zubarev and Novikov 0, E9[ 
[2(j| , |2l| and earlier papers by [^) and Cohen Q| , and cross from nonequilibrium distribution functions 
f s (x s ;t) to irreducible distribution ones G s (x s ;t), which may be introduced by equalities presented in 
^ P . In our case with some modifications we obtain 

f 1 (x 1 ;t) = G 1 (x 1 ;t), (3.2) 
f 2 (xi,x 2 ;t) = G 2 (x!,x 2 ;t) + g 2 {r\, r 2 ;t)Gi(xi;t)Gi(x 2 ;t), 

f 3 (xi,x 2 ,x 3 ;t) = G 3 (xi,x 2 ,x 3 ;t) + } j G 2 (xi,x 2 ;t)Gi (x 3 ; t) + 

p 

53(^1, r 2 , r 3 ; t)G 1 (x 1 ;t)G 1 (x 2 ;t)G 1 (x 3 ; t), 



Here, the position-dependent quasiequilibrium distribution functions g 2 (ri, r 2 ; t), g 3 {ri, r 2 , r 3 ; t), g s (r s ; t) 
are defined by relations of paper The modification of group expansions (3.2) consists in that a con- 



siderable part of the space time correlations is accumulated in quasiequilibrium functions g s (r s ;t). If all 
g s (r s ;t) = 1 for s = 2, 3, . . . these group expansions coincide with those of paper H. As far as each line 
in fl3.2| ) brings exactly for new functions G s (r s ;t), s = 1,2,3,..., the corresponding equations can be 
solved with respect to irreducible distribution functions and we may write the following: 

Gi(n;t) =fi(xi;t), (3.3) 
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G 2 (x 1 ,x 2 ;t) = f 2 (x 1 ,x 2 ;t) - g2(ri,r 2 ',t)fi(xi;t)fi(x2;i), 

G 3 (x 1 ,x 2 ,x 3 ;t) = f 3 (xi,x 2 ,x 3 ;t) f 2 (x!, x 2 ;t)f 1 (xs;t) 

p 

h 3 (r l7 r 2 , r 3 ;t)f 1 (x 1 ;t)f 1 (x 2 ;t)f 1 (x 3 ;t). 



In (3J2) and (3J3) the symbol denotes the sum of all different permutations of coordinates for three 

p 

and more particles, 

h 3 (r 1: r 2 ,r 3 ;t) = g 3 (r 1 , r 2 , r 3 ; t) - g 2 (ri, r 2 ; t) - g 2 (r u r 3 ; t) - g 2 (r 2 , r 3 ; t) (3.4) 
= h' 3 (r 1 ,r 2 ,r 3 ;t)-2, 

where h' 3 (ri,r 2 ,r 3 ;t) is the three-particle quasiequilibrium correlation function. Now let us write the 
BBGKY hierarchy [jl], |) for irreducible distribution functions G s (x s ;t), namely, the first two equations, 



-+LL(1) ) G 1 (x 1 ;t) + 



(3.5) 



J dx 2 iL(l,2)g 2 (r 1 ,r 2 ;t)G 1 (x 1 ;t)G 1 {x 2 ;t) + J dx 2 iL(l,2)G 2 (x 1 ,x 2 ;t) 



0. 



Differentiating the relation for G 2 {x\,x 2 \t) in (3.3) with respect to time and using the second equation 
from the BBGKY hierarchy for the function f 2 (xi, x 2 ; t), we can get for the pair irreducible distribution 
function G 2 {x\ ,x 2 ;t) the equation, which reads 

e) G 2 (xi,x 2 ;t) = - ( +iL 2 ) g 2 (ri,r 2 ;t)Gi{xi;t)Gi(x 2 ;t) - (3.6) 



dt 



\L 2 



f Ax 3 {ii(l,3) + i£(2,3)} {G 3 (x 1 ,x 2 ,x 3 ;t)+J2G2(x 1 ,x 2 ;t)G 1 (x 3 ;t) + 
J P 

g-i{r i , r 2 , r 3 ; t)Gi (xn t)Gi (x 2 ; t)Gi (x 3 ; t) j . 

In a similar way, we can obtain other equations for the three-particle irreducible function G3(a;i, x 2l x 3 ; t) 
and higher G s (x s ]t) ones. One remembers now, that appearing the quasiequilibrium distribution func- 
tions g 2 (ri,r 2 ;t), g 3 (ri,r 2 ,r 3 ;t), g s (r s ;t) in the hierarchy is closely connected with the fact that the 
boundary conditions for solutions of the Liouville equation take into consideration both the nonequilib- 
rium character of one-particle distribution function and local conservation laws, that corresponds to a 
consistent description of kinetics and hydrodynamics of the system jl], |22j . Since in the present paper 
we analyze two first equations ( |3.5| ) and ([T^) only, we will not write down next ones. It is important 
to note, that if we put formally g s (r s ; t) = 1 for all s = 2, 3, . . . in (3^) and (3J3), we come to the first 
two equations of the BBGKY hierarchy for irreducible distribution functions Gi{x\]t) and G 2 {x\, x 2 ; t), 
which w ere obtained in the paper [p| by D. N.Zu barev and M.Yu.Novikov. The first term in the right-hand 
side of (3J3) is a peculiarity of (|3.5| ) and (3^) equations system. This is the term with time derivative 
of the pair quasiequilibrium distribution function g 2 (r%, r 2 ;t). As it was shown in jl], |2^, ^3|, the pair 
quasiequilibrium distribution function is a functional of local values of the temperature (3{r; t) and the 
mean particle density n(r;t). Thus, time derivatives of g 2 (ri, r 2 \(3(t), n(t)) will be conformed to (3(r;t) 
and n(r;t). These quantities, in its turn, according to the self-consistency conditions p3f , will be ex- 
pressed via the average energy value ^f'(r)^ in an accompanying reference frame and via (fi(r))*, which 
constitute a basis of the hydrodynamical description of nonequilibrium state of the system. 

If we neglect the term in the right-hand side of equation (3.6), which takes into account ternary 
correlations between particles, we shall obtain from ( |3.6| ) the equation in the so-called "pair" collisions 
approximation with the following form: 



9 ' T 



: | G 2 (x 1 ,x 2 ;t) = - ( — + \L 2 ) g 2 (r 1 ,r 2 ;t)Gi{x 1 ;t)G 1 (x 2 ;t). 



The formal solution to this equation reads 



G 2 (xi,x 2 ;t) 



, Je + \L 2 ){t-t') 



df e 



dt' 



g2(ri,r 2 ;t')G 1 (x 1 ;t')G 1 {x 2 ;t'). 



(3.7) 



(3.8) 
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Inserting this solution into the first equation ( p.5[ ) we obtain the following kinetic equation for one-particle 
distribution function f\{xi]t) = Gi{xx)t): 



d_ 

at 



+ iL(l))f l (x l ',t) 



dx 2 iL(l,2)g 2 (r 1 ,r 2 ;t)fi(xi;t)f 1 (x 2 ;t) 



(3.9) 



J dx 2 i£(l,2) J dt' e 



(e + iL 2 )(t-0 
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52 (r i , r 2 ; t')fi (n ; (x 2 ; 



where the first term in the right-hand side is a generalization of the Vlasov mean field and corresponds 
to kinetic mean field theory - KMFT [^7], ^8) . The kinetic equation ( |3.S| ) is completely equivalent to that 
of which is obtained in the "pair" collisions approximation. 

Let us consider now the set of equations (3J5) and (|3.6|) in such an approximation that the three- 
particle ir red ucible distribution function G 3 (x±, x 2 , x 3) t) and h 3 (ri, r 2) r 3 ; t) are neglected in the second 
equation ( |3.6| ). This is in the spirit of the polarization approximation introduced for obtaini ng t he Lenard- 
Balescu kinetic equation for a homogeneous Coulomb plasma |l5|| . Taking into account ( 3.4) as well as 
the relation G^xi, x 2 , x 3 ; t) = and similarly h 3 (r\, r 2 , r 3 ; t) = 0, we rewrite the equation (3J3) in the 
form: 







\L 2 + e) G 2 (xi,x 2 ;t) 



d 



+ iL 2 g 2 (r 1 ,r 2 ;t)Gi{x 1 ;t)Gi(x 2 ;t) 



(3.10) 



dt ' J ' v ' \dt 

dx 3 iL(l,3)^G 2 (xi,x 2 ;t)Gi(x 3 ;t) + G 2 (x 2 , x 3 ;t)G\{xr,t)\ 

dx 3 iL(2,3)^G 2 (xi, x 2 ;t)Gi(x 3 ;t) + G 2 (xi,x 3 ;t)Gi(x 2 ;t)\ 

dar 3 {iL(l, 3) + iL(2, 3)}{g 2 (r 1 ,r 2 ;t) + g 2 (n,r 3 ;t) + g 2 (r 2 ,r 3 ; t^G^xv, f)Gi(x 2 ; t)Gi(ar 3 ; t). 

Next, let us introduce the operator, which can be obtained by variation of the Vlasov collision integral 
near a nonequilibrium distribution Gi(xx;t): 



S[ I dx 3 iL{l,3)G 1 (x 3 ;t)G 1 {x 1 ;t) 
dx 3 iL{l,3)G 1 (x 3 ;t)5G 1 (x 1 ;t) 



(3.11) 



C(xi; t)SGx(xi;t). 



Then one represents the equation (3.10) with the help of operator C{x\\t) in the form 
d 



dt 



+ IL 2 + C{xx,x 2 ;t) +e G 2 (xi,x 2 ;t) 



d 



(3.12) 



dt 



+ iL 2 + C(xi,x 2 ;t) g 2 (ri,r 2 ;t)Gi(xi;t)Gx(x 2 ;t), 



from which the formal solution for the irreducible two-particle distribution function reads 



G 2 (x 1 ,x 2 ;t) = - I dt' e^'-^Ufat') x (3.13) 

— OO 

(Jp + iL 2 + C{xi,x 2 ;t')j g 2 {ri,r 2 ,t')Gi [xi;t')Gi{x 2 ;t'), 
and U(t, t') is the evolution operator, 

U(t,t') =exp + |~/ d *" ^ 2 +£Oi,Z2;t"))J , C(x 1 ,x 2 ;t)=£{x 1 ;t)+£(x 2 ;t). (3.14) 

As a result, we obtain the expression for the irreducible quasiequilibrium two-particle distribution function 
G 2 (x\, x 2 \ t) in the generalized polarization approximation. Inserting this expression (3.13) into the first 
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equation of the chain ( |3.5| ) yields: 

(^ i +iL(l) S jG 1 (x 1 ;t)+ J dx 2 iL(l,2)g 2 (r 1 ,r 2 ;t)G 1 (x 1 ;t)G 1 (x 2 ;t) = (3.15) 
t 

J dx 2 J dt' e^'-'hLil, 2)U(t, f ) (JL + \L 2 + C( Xl ,x 2 ; i')) g2(r u r 2 ;t')G 1 (x 1 ;t')G 1 (x 2 ;t'). 

— oo 

This is the kinetic equation for the nonequilibrium oneparticle distribution function with the non- 
Markovian collision integral in the generalized polarization approximation. It is necessary to note, that 
the presence of the Vlasov's operator C(xi,x 2 ;t) in the collision integral ( 3.15 ) indicates about taking 
into consideration collective effects. An analysis of the collision integral (3.15) in the general case is 
rather a complicated problem. But it is obvious, that the collision integral in ( 3.15| ), or the expression 
for G 2 (xi,x 2 ;t) in ( 3.13| ) may be much simplified fo r eve ry physical model of a particle system, or for 
each nonequilibrium state of the collision integral in (3.15). To show this, we shall consider two concrete 
cases: a hard spheres model and a Coulomb plasma. 



4 Hard spheres models in the polarization approximation 



In this section we shall perform the investigation of kinetic processes for hard spheres model in ap- 
proximations, which are higher than "pair" collisions one. We take into account the character of model 
parameters and the results of the previous section of this article and papers fl4 , 15 29). This investigation 
is convenient to carry out on the basis of equation chain (3.5), (3.6) at formal substitution of a potential 
part of the Liouville operator iL(l, 2) ||, |2) by the Enskog collision operator T(l, 2) jj], 0. In this case, 
equations ( |3.5| ) and ( |3.6| ) have the form 

(4.1) 



-+LL(l))G 1 (x 1 ;i) + 



dx 2 T(l, 2)g 2 {r u r 2 ;t)G 1 {x 1 ;t)G 1 {x 2 -t) + / dx 2 f (1, 2)G 2 ( Xl ,x 2 ;t) = 



^+iL§ + f(l,2)+e) CVo(. ,-,..,■;:/) 



(4.2) 



^+ii§ + f(l,2)^ fl2 (ri,r 2 ;t)(3 ! i(a;i;t)Gi(a: a ;t)- J dx 3 {f (1, 3) + f (2, 3)} 
^G 3 {xi 1 x 2 ,x 3 ;t) + ^ G 2 (x 1 , x 2 \ t)Gi(x 3 ; t) + g 3 {r 1 ,r 2 ,r 3 ;t)Gi(x 1 ;t)Gi{x 2 -t)G 1 (x 3 ;t)Y 



Further, we will consider the same approximations concerning to equat ion ( 4.2), in which G 3 (x\, x 2 ,x 3 ; t) 
and h 3 (ri,r 2 ,r 3 ;t) are neglected. Then, if we introduce similarly to (3.11) the Boltzmann-Enskog colli- 
sion operator C(xi;t), the equation (4.2) could be rewritten in the next form: 



6 / dx 3 f(l,3)G 1 {x 1 ;t)G 1 (x 3 ;t) = C(x 1 ;t)5G 1 (x 1 ;t), 



(4.3) 



Of 



+ iL° + f (1, 2) + C{x u x 2 -t) + e G 2 {x u x 2 ;t) 



d 



dt 



+ iL° 2 + T(l,2) + C{x 1 ,x 2 ;t))g 2 {r 1 ,r 2 ;t)G 1 (x 1 ;t)G 1 (x 2 ;t), 



Hence it appears, that the formal solution to G 2 {x\, x 2 ; t) reads 



(4.4) 



G 2 (xi,x 2 ;t) 



dt' e^'-^t/hsM') x 



(4.5) 



— + ill + f (1, 2) + C( Xl ,x 2 ;t') \ g 2 {n,r 2 -J)G x { Xl] t')G x {x 2 -J), 
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where Uh s (t,t') is the evolution operator for the system of hard spheres: 



U hs (t,t') =exp + { - / dt" \L° 2 + T{l,2) + C{x u x2;t") 



(4.6) 



C(xi,x 2 ;t) = C(xi;t) + C(x 2 ;t). 



Now let us put (4.5) into the first equation (4.1). Then the resulting equation looks: 
d 



dt 



LL(1) Gi{xi-t) = / dx 2 T(l,2)g 2 (r 1 ,r 2 ;t)G 1 (x 1 ;t)G 1 (x 2] t) - / dx 2 T(l,2) 



(4.7) 



dt' e' 



^'-^U h3 (t, t') | A + iL" + f (1, 2) + C( Xl ,x 2 ; t') j ff2 (n, r 2 ; i')Gi (hjOGi (a*; i 7 )- 



This equation can be called as a generalized kinetic equation for the nonequilibrium one-particle distri- 
bution function of hard spheres with the non-Markovian collision integral in the generalized polarization 
approximation. The first term in the right-hand side of this equation is the collision integral from the re- 
vised Enskog theory |30|. Neglecting time retardation effects and assuming that the operator C(x\, x 2 ; i) 
does not depend on time when 



Gi(<Bi;t) = /o(p) 



2nkT 



3/2 



exp 



P 



2mkT 



is the local equilibrium Maxwell distribution function, the next term can be rewritten in a simplified 
form, 

* t 
I R (xi;t) = - J dt'e^'-^oOn^OGiOn;*')- J df e e ( t '-^R 1 (x 1 ;t,t')G 1 (x 1 ;f), (4.8) 



where 



R (x i; t, = f dx 2 f (1, 2)e{( t '- t X iL " +f (i.2)+c(*i,* a ))}[ iL o +C (x 1 ,x 2 )] g 2 (n, r 2 ; t^G^t'), (4.9) 



R 1 (x 1 ;t,t') = [ dx 2 f(l,2)e{( t '- t K i ^+t(i^)+c(x 1 ,, 2) )} f(lj2)52(ri;r2 .^ )Gi(x2 .^ ); (4JQ) 



Rx{xx\t,t') is the generalized ring operator. The kinetic equation Q4.7| ) together with (|4.9|), ( [4.10[ ) is the 
generalization of the kinetic equation for a system of hard spheres, which has been obtained by Bogolubov 
in [ p"5"[ p9| . It coincides with that, when the quasiequilibrium pair distribution function of hard spheres 
is set formally to be unity. 



5 Coulomb plasma in the polarization approximation 

Here we shall study an electron gas, which is contained into a homogeneous positively charged equilibrat- 
ing background. This background can be created, for example, by hard motionless ions. Then, electrons 
interact according to the Coulomb law: 



e 1 



\r\-r 2 \ |ri2| 

the Fourier transform of which exists in the form of a real function #(|fe|): 

f__t _dk_ w ik . ri2 _ w 

ri 2 ~ J (2tt)3 fc2 e . *W - k 2 . {y- 1 ) 

here k is a wavevector, e is the electron charge. Let us consider equation chain ( |3.5[ ), ( |3 . 10| ) in the homoge- 
neous case, when G\{xt\t) = Gi(p 1 ;t) and pair distribution functions depend on \ri 2 \. Following the Bo- 
golubov method Jl5l, [III , we shall suppose that the one-particle distribution function G\ {p 1 ; t) is calculated 
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in the "zeroth" order on interaction constant q, pair distribution functions G2(ri 2 ,p 1 ,p 2 ] t) and <?2( r i2; t) 
as the first order on q, and G3(xi,x 2 ,X3;t), gz(r\,r2,rz;t) ~ q 2 , where q = ^-0, r<j = ^/ (9/47re 2 n is 
the Debye radius, n = N/V, = k^T, fee is the Boltzmann constant, T is thermodynamic temperature. 
Therefore, to obtain an equation for G 2 (ri 2 , Pi, p 2 ; t) in the first approximation on interaction constant 
q without time retardment effects, it is necessary to retain all integral terms, but omit all the others. In 
this case, using the Fourier transform with respect to spatial coordinates for a homogeneous Coulomb 
electron gas, the set of equations (3.5), (3.1C) yields: 

^G x {p x] t) = J dk&p 2 i${\k\)g 2 {k^^ 

^Gi(Pi;t) = ^-G?i(Pi;t) J dk k-P(\k\)^mg 2 (k;t) + A J dk k$(\k\)ZmG 2 (k, Pl ;t) (5.2) 

and equation for G 2 (k,p 1: p 2 ;t): 
' d_ 

at 



or 



(5.3) 



ik$(\k\) 
ik$(\k\) 



-Gt(pi,t) J dp 3 G 2 (k,p 2 ,p 3 ;t) - ^-Gi(p 2 ;t) J dp 3 G 2 (k,p v p 3 ; 



_d 
dp 

— G 1 (p 1 ;t) 52 (-fe;t)Gi(p 2 ;i)-^-G 1 (p 2 ;t) ff2 (fe;t)G 1 (p 1 ;t) 



t) 



e — > +0, and G 2 (k 7 p 1 ;t) = J dp 2 G 2 (k,p 1 ,p 2 ;t); ^smg 2 (k;t), 3m G 2 (x\, x 2 ; t) are imaginary parts of 
the corresponding distribution functions. The following properties should be noted: 



G 2 (-k,p 1 ,p 2 ;t) 
92(-k;t) 



G* 2 (k,p 1 ,p 2 ;t), 
g*(k;t). 



The solution to (5.3), neglecting time retardment effects, reads: 

k<P(\k\) 



G 2 (fc,Pi,p 2 ;t) 



P| , i^Gi(Pi;t)G 2 (-fe,p 2 ;t)-— G 1 (p 2 ;i)G 2 (fe,p 1 ;t)|+ (5.4) 
m 

p! ' J — G 1 (p 1 ;t) 52 (-fe;t)G 1 (p 2 ;t) - — GifaMk^G^t) 

k . til _ i0 I <9pi <9p 2 



It should be noticed also that equation (5.2) contains an imaginary part of the pair irreducible nonequi- 
librium distribution function, to be integrated with respect to momentum of the second particle. Now 
one integrates equation (5.4) over all values of momentum p 2 and defines in such a way some function 
G 2 (fe,p 1 ;t): 



1 



_d 
dp 



d P2 ^-G,(p,:l) 1 G L ,(fc.p,:/l 

k ■ £ - i£ - iO °Pi 

m 

-G x (p i; t) f dp 2 fc 5J fc|) G 2 (-k,p 2 ;t) 
m 

fc$(|k|) 



(5.5) 



... [ 9 9 

dp2 fc P12 r |^- G i(Pi; < )52(-fc;0Gi(p 2 ;i) - ^G f i(p 2 ;% 2 (fe;t)Gi(p 1 ;t) 

m 

Further, we should exclude from ( |5.5| ) the term with G 2 (— fe,p 2 ;f). To do this, we follow Lenard [[33], Q 
and integrate the equation ( |5.5| ) over momentum component p^, which is perpendicular to wavevector 
k. Resulting expression then reads: 



8 f k<P(\k\) 

[l + $(\k\)X(k,px;t)]G 2 (k, Pl ;t) = — G 1 {p 1 ;t) / dp 2 " » G 2 (-k,p 2 ;t) 

°Pi J k-'—-iO 

m 

J Am k$(\k\) / d 



(5.6) 



dp 2 



k ■ 



, G 1 (pi;i) 52 (-fe;i)Gi(p 2 ;t)--^-G 1 (p 2 ;i) 52 (fe;t)Gi(p 1 ;t)^ 
W2 iO V d Pi d P2 
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Here the following conventional designations have been introduced: 
X(k,pi,t) = / dp 2 



k .Pll. i0 dp 2 
m 

G 1 (p 1 ;t) = J dp lx G x (pi;t), 
G 2 (k,p 1 ;t) = J dp 1± G 2 {k,p 1 ;t). 



k Q fy ' k f) ' k 

Gi(p 2 ;t), pi = — — , p 2 = — , k = \k\ 



(5.7) 



Now we multiply both equations ( pT5| ) and (5.6) on ■^■Gi(pi;t) and ■^■G±(p 1 ;t), respectively, and 
subtract them: 



(l + #(|fc|)X(fc,Pi;t) 
$(\k\)X(k, Pl ;t)g 2 (k;t) 



G 2 (k, Pl ;t) — Gx (pi;t)-G 2 (k,pi;t)—G 1 (p 1 ; t) 
dpi dpi 

Gi (pi ; t) — G x (p x ; t) - G x ( Pl ;t)—Gi ( P i ; t) 
dpi dpi 



(5.8) 



If we extract imaginary part of this equation, one can find unknown quantity ( 3mG 2 (k,p 1 ; t), provided 
9fmG 3 (*,Pi;t) = " 



^-Gi{pi;t) 3raG 2 (fc, Pl ;i) = 
dpi 



<g(|fc|)3m [X(fc,Pi;t)g 2 (fc;t)] 
|l + <2>(|fc|)X(fc, Pl ;i)| 2 

Gi (pi ; t) — G x ( Pl ;t)-Gi(p 1 ;t)—Gi (pi ; t) 



(5.9) 



Since 3mX(fe,pi;i) = —n-^Gi(pi\t) |3j| putting the expression for Shn G 2 {k, p x ; t) into the equa- 
tion (5.2) gives the generalized Bogolubov-Lenard-Balescu kinetic equation for an electron gas in an 
equilibrating background 



Q-Gi{Pi,t) = —Gi( Pl ;t) J dk k<P(\k\)Zmg 2 (k;t) 
) f 

— J dp 2 Q(pi,p 2 ;t) 



(5.10) 



_d_ 

dpi 



d d 
dpi dp 2 



Gi( Pl ;t)Gi(p 2 ;t), 



where Q(Pi,p 2 ',t) is a second rank tensor 



Q(Pi,P2\t) = / dk — 
J 1 



\$(\k\)rk- k 



$(\k\)X(k,pi;t)[ 



■$smg 2 (k;t)6(k- (p x -p 2 )), 



(5.11) 



which coincides with Q(p x , p 2 ) [|32| at ?img 2 (k; t) = 1. In this case, the kinetic equation ( 5.10 ) transforms 
into the well-known Lena rd-B alescu equation |32[ [33| |lj . Evidently, the generalized Bogolubov-Lenard- 
Balescu kinetic equation ( 5.10 ) climes the description of a dense electron gas, since in the both generalized 
mean field and generalized Bogolubov-Lenard-Balescu collision integrals, manyparticle correlations are 
treated b y ima ginary part of g 2 (k; £)• Nevertheless, the problem of divergence in the collision integral of 
equation (5.10) at small distances (k — ► 00) still remains. There are papers, where divergence of collision 
integrals is avoided with the help of special choice of differential cross section (quantum systems |35| ) , or 
via combination of simpler collision integrals (classical systems jn|). These generalization for collision 
integrals are attractive by theirs simplicity and are usable for ideal plasma. But, contrary to the obtained 
by us Bogolubov-Lenard-Balescu kinetic equation, they do not work for nonideal plasma. In accordance 
to the proposed structure of collision integral 



l total 



-feoltz 



andau 



cnard— Balcscu 



the influence of particles interaction into plasma energy will be defined by correlation function 52(f)- Its 
asymptotic is 



lim 



exp 



e a eb 
rk B T 



e a eb 



Boltzmann 



rk^T rkftT 

Landau 1 



e a e b c -r/r D 



d — Balcscu 



(!), where r D denotes the Debye radius. 
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At other combinations one arrives at false expressions for thermodynamic functions |13| ]. Dynamical 
screening, appearing in the obtained by us generalized Bogolubov-Lenard-Balescu collision integral, is 
free of these discrepancies. Generally speaking, the problem of divergency could be solved within a frame 
of charged hard spheres model, combining the results of this section and previous one. But this step 
constitutes an intricate and complicated problem and needs a separate consideration. 

Evidently, an investigation of the obtained kinetic equation is important in view of its solutions and 
studying transport coefficients and time correlation functions for model systems. 

6 Conclusion 

In view of dense systems study, where the consideration of spatially interparticle correlations is important, 



the BBGKY hierarchy (p.5|), (3.6) with the modified boundary conditions and group expansions has a 



quite good perspective. The kinetic equation (4.7) - (4.10) is a generalization of the Bogolubov one 
|15| , |29fl for a system of hard spheres. M.Ernst and J.Dorfman |3l| had investigated collective modes in 
a nonhomogeneous gas and showed, that the solution of a dispersion equation for hydrodynamic modes 
leads to the nonanalytic frequency dependence on wavevector. This is connected with the fact, that 
the ring operator for nonhomogeneous systems at small wavenumbers has a term proportional to Vk. 
Similar investigations of collective modes and time correlation functions in the hydrodynamic region 
have been carried out by Bogolubov jlj| |||, Nevertheless, it is necessary to carry out analogous 
investigations of hydrodynamic collective modes and time correlation functions on the basis of kinetic 



equation fl4.7|), taking into account ( |4.8| ) - fl4.10| ), where some part of space correlations is considered in 
the pair quasiequilibrium function 52( r i, r 2] t). Obviously, these results may appear to be good for very 
dense gases, which could be described by a hard spheres model. An important factor is that in the kinetic 



equation (4.7) - ( 4.10| ) as well as in the generalized Bogolubov-Lenard-Balescu one, collective effects are 



taking into account both via Vlasov's mean field and pair quasiequilibrium correlation function, which is 
a functional of nonequilibrium values of temperature and chemical potential. 

Transferring the obtained results on quantum systems is not obvious. Such a procedure is rather 
complicated and needs additional investigations. Nevertheless, some steps in this way have been done 
already by our colleagues |36|, |37j . 
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